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^ I Abstract 

I— I ■ The lilypond model on a point process in d-space is a growth-maximal system of non- 
overlapping balls centred at the points. We establish central limit theorems for the 

I total volume and the number of components of the lilypond model on a sequence of 

• Poisson or binomial point processes on expanding windows. For the lilypond model 

d . over a homogeneous Poisson process, we give subexponentially decaying tail bounds 

^ I for the size of the cluster at the origin. Finally, we consider the enhanced Poisson 
lilypond model where all the balls are enlarged by a fixed amount (the enhancement 

T— I ■ parameter), and show that for d > 1 the critical value of this parameter, above which 

I the enhanced model percolates, is strictly positive. 
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2 : 1 Introduction 



X 



Suppose (f is a. locally finite set of points of cardinality at least 2 in the space M"^. The 
lilypond model based on is the system of balls (or grains) {Bp(^x){x) : x G (here 
5^ I Br{x) := {y E M.'^ : \y — x\ < r} and | ■ | is Euclidean norm) with the following two 

properties: 

• The hard-core property: p{x) + p{y) < \x — y\ for all different x,y & (p. 

• The smaller grain-neighbour property: for each x G there is at least one y G (p\{x} 
such that p{x) + p{y) = \y — x\ (in which case the points x and y are called grain- 
neighbours) and p{y) < p{x). In this case we call y a smaller grain-neighbour of 

X. 
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In the sequel we shall write p(a;, ip) to denote the dependence of the radii on both ip and 
X G v^. Heveling and Last [5] established existence and uniqueness of the model for all 
such if (in fact in greater generality). Define the union set 

Z(y,):= Ui?p(x,^)(a;)- (1-1) 

In the case of finite if, the lilypond model may be constructed as follows. All points of 
if start growing at the same time and at the same rate. Any given ball ceases its growth 
as soon as it reaches any other ball. When has just a single element x, we define 
p{x, ip) := +00. 

The model is of interest since it is a growth-maximal hard-core model. No single ball 
can grow further without overlapping another ball. It also has a maximin property: if 
ip is finite with n elements and the radii {p(x, y?) : x G are listed in ascending order 
as pi,p2, ■ ■ ■ Pn, and if {p'{x) : x G 9?} is any other system of radii satisfying the hard- 
core property, similarly listed in ascending order as p[, ■ ■ ■ , p'n then (pi, . . . exceeds 
{p[, . . . , p^) in the lexicographic ordering on M". 

In this paper we consider the lilypond model on random ip. The lilypond model was 
introduced by Haggstrom and Meester [8] for the case where is a stationary Poisson 
process $ in of intensity one; we call this the Poisson lilypond model and set \1/ = 
{(x,p(x, $)) : X G $}. They proved that the union set Z := does not percolate, 

i.e. does not have an unbounded connected component. Interestingly, there does exist a 
stationary percolating hard-core system of (non-lilypond) grains on $, at least in high 
dimensions; see [2]. 

Apart from the one-dimensional case (see Daley, Mallows and Shepp [1]), only a few 
further probabilistic properties of \1/ are known. Daley, Stoyan and Stoyan [5] give bounds 
for the volume fraction 

pz := P(o ez)= EV,(Z n [0, l]'^) (1.2) 

of Z, where denotes the origin in and Vd is Lebesgue measure (volume) on M'^. The 
latter paper also has some numerical results on the typical radius po of The distribution 
of Po is that of p(x, $) for a "randomly picked" x G Because $ is Poisson, it is well- 
known that the distribution of po is that of p(0, $°) where $° := $ U {0}. Because of the 
hard-core property of ^ we clearly have 

Pz = hdE[ptl (1-3) 

where hd := Vrf(i?i(0)) is the volume of the unit ball in W^. 

The contributions of the present paper fall into three categories: tail bounds, central 
limit theorems, and non-percolation under positive enhancement. These may be viewed 
as extending the percolation theory of the Poisson lilypond model beyond the basic fact 
that Z does not percolate. We now give an overview of our results, which are proved 
using notions of stabilization developed in Sections [2] and |3l 

Given ip> as above, and given x G v^, define C"(x, y?) to be the connected component 
of Z[lp) containing x. We study the 'cluster at the origin' C"(0, $°), which amounts to 
studying C"(x, $) for a randomly picked point x G Since Z does not percolate, we 
know that C"(0, $°) is almost surely bounded. In Section H] we shall give tail bounds on 
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its size. We consider three different measures of the 'size' of C"(0, $'^), namely volume, 
metric diameter, and the number of constituent grains. In each case, we shall give upper 
tail bounds showing that the probability of the size of a cluster exceeding t decays expo- 
nentially in a power of t. We also give lower bounds of the same form but with different 
exponents; it remains open to establish the 'correct' exponent (if any) for the tail decay 
of the size of the cluster at the origin. 

Our central limit theorems are stated with reference to a sequence of expanding win- 
dows in W^. Let C M'^ be a compact convex set containing an open neighborhood of the 
origin with VdiW) = 1, and set Wn := n^l'^W , where for A C M'^ and x G M*^, a G M+ we 
write X -|- aA for {x + ay : ?/ G A}. For A C M*^, let denote the number of connected 
components of A. In Section [5l we shall derive central limit theorems for (among other 
things) Vd(Z($„)), and fi;(Z($„)), where we set := $nPKi. We shall also estabhsh the 
corresponding de-Poissonized central limit theorems, where instead of one considers 
a point process Xn consisting of n independent uniformly distributed points in n^^'^W. 
Central limit theorems such as these could be of use in establishing confidence intervals 
for quantities such as the volume fraction, based on simulations. 

In Section E] we consider the enhanced lilypond model Z^, defined for 5 > to be the 
^-neighbourhood of Z in M*^. We shall show that for d > 2, there is a strictly positive 
critical value of 6 such that the enhancement Z^ percolates almost surely if 6 is above 
this critical value, and does not percolate if 6 is below the critical value. This gives us 
a one-parameter family of random sets in the continuum exhibiting a non-trivial phase 
transition. 

As a final remark here, we compare the Poisson lilypond model to the random sequen- 
tial adsorption (RSA) process with infinite input, in which unit balls arrive at locations 
given by a homogeneous space-time Poisson process starting at time zero, each ball being 
irreversibly accepted if it does not overlap any previously accepted ball. Continuing to 
time infinity, one ends up with a maximal system of balls satisfying the hard-core con- 
straint, as with the lilypond model; see Penrose [12], and Schreiber, Penrose and Yukich 
[IB] for formal definitions, existence and limit theorems for RSA. Noteworthy differences 
are that for RSA, the radii are all the same, the point process of accepted ball centres is 
not a spatial Poisson process, and the balls almost surely do not touch. 

Notation: We use c, c', c" and so on to denote positive finite constants whose values 
are unimportant and may change from line to line. On the other hand, we denote by Ci, C2 
and so on, constants whose values (though still not very important) are fixed and which 
may reappear in other parts of the paper. For nonempty A C M'' we write diam(A) for 
^'^P{x,y)<£AxA k " y\- We write card(A) for the number of elements of A (possibly oo). 

2 Stabilization 

In this section, we establish that there is an almost surely finite random variable R : = 
-R($) such that the radius p{0, $ U {0}) is unaffected by modifications to the point set $ 
outside the ball B^{0). This is known as a radius of stabilization for p(0, $). Moreover, 
we establish tail bounds for R, i.e. bounds on the ¥{R > t), which decay exponentially 
in . These tail estimates will be crucial in all of our subsequent proofs. 

A point process is defined as a random variable taking values in the space N of all 
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locally finite subsets of equipped with the smallest cx-field A/" containing the sets {if G 
N : (p{B) = k} for all Borel B G M.'^ and all k G Nq, where f{B) denotes the number of 
elements of (f (1 B. By a (finite) descending chain in G N we mean a finite sequence 
xo,---,Xn {n > 1) of distinct points of ip for which \xi-i — Xi\ > \xi — Xi+i\ for all 
i G {l,...,n — 1}. Note that any two points of (f, considered on their own, form a 
descending chain. Let G N and x E (p. If (y9(M'^) > 2 we let D{x,ip) denote the nearest 
neighbour distance from x in ip\{x}. That is, we set D{x,ip) := min{|x — ?/| : y G ip\{x}}. 
U ip = {x}, set D{x, if) := +oo. 

For any G N and x,y G M!^ we set ip^ := if VJ {x} and y?^''^ = U {x,y}. We 
construct a closed set Siy, ip) C such that if this set is bounded, then the radius 
p{y,(py) is determined by the restriction of ip to S{y,(p) (see (12.21) below). In the trivial 
case where ip \ {y} = we define S{y, (p) := M'^. Otherwise, we define 

S{y,^):=B2Diy,^y){y)^ U 5,(x), (2.1) 

{x,r)i^A{y,ip) 

where the set A{ii, (p) G ip x (0, oo) is defined as follows. A pair {x, r) belongs to A{ip>) if 
there is a descending chain xq, . . . , a;„ in ipP such that xq = y, \xi —y\< 2D{y, (p^), Xn = x 
and r = \x — 

The next result tells us essentially that S{y, $) is a stopping set; see [1], [T7] . 

Lemma 2.1. Suppose ?/ G M'' and G N wit/i ■?/' fl S{y,ip) = (p> H S{y,(p>). Then 
S{y,i^) = S{y,ip). 

Proof. Assume (p \ {y} is nonempty (otherwise the result is trivial). Then D{y, (p^) = 
Diy.r)- 

Suppose (x, s) G A{y, ip). Then there is descending chain xq, . . . , x„ in ip^ such that 
Xo = |xi — y\ < 2D{y,ipy), x„ = x and s = |x — x„_i|. For 1 < m < we have 
Xm G S{y, p>) Hip so that Xm £ "^A- Hence, (x, s) G v4(?/, ?/'), so A{y, ip) C v4(?/, ip). 

Conversely, suppose (x, s) G A{ii, ip). Then there is descending chain xq, . . . , x„ in ip^ 
such that Xo = y, |xi — y\ < 2D{y,ipy), x„ = x and s = \x — x„_i|. 

We claim that x^ G for 1 < m < n. This is proved by induction on m; if it holds 
for 1 < m < A; — 1 then x^ G S{y, ip) (1 tp so Xk G To start the induction note that 
\xi — y\ < 2D{y, ipy) so Xi G Siy, ip) P\ ip so Xi E ip. 

By the preceding claim, (x, s) G A(?/, v?) so y4(y,'?/;) C Aiij^ip), and hence Siy^'p) C 
S{y,ip). Therefore Ad/,-?/;) = A(y,y9), so that Siy^ip) = S{y,ip). □ 

Lemma 2.2. For any ?/ G M'^ anc? ip eN, if the set S{y, ip) is bounded, then it satisfies 

p{y, ipy) = piy, {ipy n S{y, ip))Ui;), V C M'^ \ S(y, y^), G N. (2.2) 

Proof. Suppose S{y,ip) is bounded and G N with ^/^ C M"' \ 5'(v9). Since S{y,ip) = 
S{y,ipy), we can and do assume without loss of generality that y E ip. Also, by Lemma 
12.11 it suffices to prove the result in the case where ip C S{y, ip), so we assume this too. 
Write ip' for ipUip, and write p(x) for p(x, ip) and p'{x) for p(x, v?')- 
Suppose p(?/) > p'{y). Let xi be a smaller grain- neighbour of y in ip'. Assume for 
now that Xi G Then by the hard-core property of the lilypond model on ip, we have 
piy) + ^ ki ~ 1/1 = p'(l/) + p'l^^i); cind hence p(xi) < p'(^i)- 
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Next let X2 be a smaller grain-neighbour of Xi in (p. By the hard-core property of the 
lilypond model on we have p'(xi) + p'(x2) < X2I = p{xi) + p{x2), so p'{x2) < p{x2). 
Let be a smaller grain-neighbour of X2 in Assuming X3 G V9, using once more the 
hard-core property of the lilypond model on ip yields p{x^) < p\x^). 

Continuing to alternate in this way, we get a sequence of points Xi satisfying 

pi.y) > p'iy) > Pi^i) > Pi^i) > Pix2) > p'{.X2) > p'ixs) > p{xs) > ■■ 

and terminating at x„ if for some (odd) n we have x„ G V'- We see from these inequali- 
ties that the possibly terminating sequence y,Xi,X2, ■ ■ ■ consists of distinct points. Also, 
setting xq = y and pi = p{xi) and p- = p'{xi), we have \xi — Xi_i\ = p'i + p'i_i for odd i 
and \xi — = Pi + Pi-i for even and therefore from the above inequalities \xi — 
is nonincreasing (in fact, strictly decreasing) in i. 

Thus the (possibly terminating) sequence |/,Xi,X2, . . . forms a descending chain with 
~ 2/1 < 2D(?/, (py). If the sequence terminates at some point Xn = z with z & if), then 
z G S{y,ip), contradicting the assumption that i/' H S{y,ip) = 0. On the other hand, if 
the sequence (xi) does not terminate, then it forms an infinite descending chain in 93, 
contradicting the assumption that S{y, ip) is bounded. 

Thus if p{y) > p'{y) we have derived a contradiction. If piy) < p'iy) we argue similarly, 
this time starting with xi a smaller grain-neighbour of y in ip>. Again we end up with a 
contradiction. Thus we must have p{y) = p'{y). □ 

Given G N, define 

R{^) := inf {r > : ^(0, ^p) C 5^(0)}, (2.3) 
with the convention inf(0) := +00. 

Lemma 2.3. The function i? : N — ^ [0, 00] is Borel-measurable. 

Proof. Let t > 0. It suffices to prove R~^{[0,t\) is measurable, i.e. in A/". For m G N, 
let = G N : (^(5^(0)) = m}. Then by Lemma [SHI 

R-\[0, t]) = U^=o{^ e : R{v n BtiO)) < t}. 

However, it is not hard to see that for all m G N, the function g : 5^(0)™ — )■ {0, 1} given 
by 

'^{R({xi,...,xm})<t} a xi,...,Xm are distinct 
otherwise 



g{xi, ...,Xm)--- 

is measurable, and hence {ip G : R{{p fl -Bt(O)) < t} G Af. The result follows. □ 



By Lemma [2.21 R{p>) is a radius of stabilization for p(0, p>^), i.e. p(0, p>^) is unaffected 
by modifications to p> outside Bpi(^^^{0). Also, -R($) is known to be almost surely finite 
[TU| 13]. The next result provides another proof of this last fact, and more importantly 
shows that the tail of the distribution of -R($) decays sub-exponentially. 

Lemma 2.4. There is a constant Ci > such that for any x G Mf^, 

F{R{^) > t) < ci exp(-q 4'^/('^+i)), t > 0. (2.4) 
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Proof: Let t > 1. The first step is to bound the distance D := D(0,$'^) from the 
origin to its nearest neighbour in $ by a suitable power of t. For the duration of this 
proof, define 

€ := l/{d +1); K:= max(2, (26de')^/('^+^)). 

Given t > 1, set u := K'^f^ and ^ := \Kt^~'^\. Let E be the event that there is a 
descending chain 0,a:i, . . . ,a;f in $° such that K\xi\ < (i.e., |xi| < u). Then we assert 
that the event inclusion 

{R{^)>t} (l{2KD>f}UE (2.5) 

holds. To see this, assume on the contrary that 2KD < and that E does not occur. 
Then any descending chain starting at and with its first link of Euclidean length < 
2D < K'^f^ would have at most Kt^~'^ — 1 links, so would end at a point of Euclidean 
norm at most u{Kt^^^ — 1). This would imply that -R($) < uKt^~^ = t. Hence (12. 5p 
holds. 

For any set A, let {A)i denote the set of ^-tuples (ai, . . . , a^) G such that ai, . . . , 
are distinct. Then, using the ith order Palm-Mecke formula for the Poisson process (see 
e.g. Theorem 1.6 of fU]) similarly to Subsection 3.2 of |3], we have that 

P(^) < E 1{m > > |a;2-a;i| > ... > Ix^-x^-il} 

{xi,...,Xf)e{'S>)i 

/.../l{.>|x.|>|x,-x,|>...>|x,-x,_,|}<i.,...<ix, 

> • • • > \yi\}dxidy2. . .dyi 



> 


\xi\ 


> \. 


> 




> 1 



{bdU 



(2.6) 

Hence by Robbins' bound associated with Stirling's formula (see e.g. [6]), and the defini- 
tion of u and £, and the fact that K > 2, 



^(^) ^ £^(27r)i/2 ^ \Kd{Kt^-e - 1)7 ^ \Kd+H^-\ 
By definition, ed = 1 — e, and K'^'^^ > 2bde^ so that 

P(^) < e-^ < exp(-(ir/2)t'^/(^+^)). 
Returning to (12.51) and using the definition of e again, we see that 

P(i?($) >t)< exp(-6d(2ir)-^t^/('^+i)) + exp(-(ir/2)t'^/('^+^)), (2.7) 
and therefore (12. 4p holds for suitably chosen Ci. □ 

Next we extend Lemma 12.41 to the family of binomial point processes 

Xn,m ■■= W'^Xi : 1 < « < m}, n, m G N. (2.8) 

where Xi,X2,... is a sequence of independent random rf- vectors uniformly distributed 
over the convex set W (which was introduced in Section [1]). 
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Lemma 2.5. There is a constant C2 > 0, dependent only on the choice of W , such that 
for all n & N with n> 4, and all x G Wn, and all m E [n/2, 3n/2], 

nni-x + Xn,m-i) >t)< C2exp(-C2-4'^/(^+i)), t > 0. (2.9) 

In proving this, and again later, we shall use the fact that there is a constant Cq > 0, 
dependent only on W, such that 

Co^r'^ < Vd{WnnBr{x)), neN,x eWn,r e (0, n^/'^diam(iy)]. (2.10) 

To see this, take m > such that -B2m(0) C W. For x G W, and < s < u, the 
convex hull of {x} U B2u{0) is contained in W, and the intersection of this with Bs{x) 
contains the intersection of Bs{x) and a cone with apex at x and subtended angle at least 
arcsin('u/diam(Vr)), and hence has volume at least c~^s'^ for some constant c independent 
of X and s. Therefore 

mi{s-'^Vd{W n Bs{x)) : X G Ty, s G (0, m]} > 0. (2.11) 

Since s~'^Vd{W fl Bs{x)) is continuous in (s, x), it is bounded away from zero on (s, x) in 
the compact set [u, diam(Ty)] x W , and therefore (12.111) still holds with the range of s 
extended to (0, diam(Ty)], and then (I2.10p follows by setting x = n~^^^y and scaling. 

Proof of Lemma \2. 51 By (12.101) . for all n > 4, x G Wn, m > n/2 (so that in particular 
m — 1 > n/4) and r G (0, n-'^/'^diam(W^)], we have 

P(D(x,Xn,™-i) > r) < (1 - Co^n-V^)™-! < (exp(-Co ^n-V^))™"! 

< exp(-(4co)-V'^), (2.12) 

and this holds trivially for r > n^/'^diam(Vl^) as well. 

The proof of Lemma 12.51 now mostly follows that of Lemma 12.41 There is a difference 
in the first term in the right side of (12. 7p . where we now need to estimate the probability 
that the nearest point to x in Xn,m-i is at a distance greater than {2K)~^t'^. For this we 
can use fl212|) . 

Instead of the event E featuring in the proof of Lemma 12. 4[ we now need to consider 
E', defined to be the event that there is a descending chain (x, Xi, . . . , x^) in Xn m-i "with 
l^^i — x\ < u. Instead of the estimate (12. 6p . setting Yi := n^^^Xi we now have 

F{E') <E Yl ^ 1^*1 1 >\y^2-y^A>■■■> \y^, - y^j} 

(il,i2,...,jf)e({l,2,...,m})^ 



/ '^{u > \xi\ > \x2 — Xi\ > ■ ■ ■ > \x£ — x^_i I }n ^dxi . . . dx£ 

< 2^ J ■ ■ ■ J l{u > |xi| > |x2 — xi| > ... > \xe — xe-i\}dxi . . . dxi (2.13) 



because m < 2n. 

With these changes, we can complete the proof by following the proof of Lemma 12. 4j 
it is easy to modify the argument to allow for the extra factor of 2^ in (I2.13P compared 
to dM]). □ 
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3 External stabilization 



In this section we introduce the notion of a fence. Loosely speaking, given an annulus in 
W^, a configuration (f has a fence if there are enough points in the annulus to guarantee 
that no lilypond grain centred inside the annulus can penetrate too far outside, and no 
lilypond grain centred outside the annulus can penetrate too far inside. Combining this 
with with the notion of a radius of stabilization R{(p) as already considered, we shall arrive 
at a stronger external stabilization radius, denoted Re^{(p), with similar tail behaviour. 
Loosely speaking, external stabilization means that changes to if beyond distance i?ex(v^) 
do not affect the grains near the origin, and changes near the origin do not affect the 
grains beyond distance Re^{ip). 

Let a; G M'' and < r < s. Let B°{x) be the open ball of radius r centred at x, and 
set dBs{x) := Bs{x) \ B°{x), the boundary of Bs{x). Let (ei, 62, 63, . . .) be an arbitrarily 
chosen sequence forming a countable dense set in dBi{0). Let k{s,r) be the smallest 
positive integer k such that there exists an increasing sequence (ji, j2, • • • ,jk) of positive 
integers, such that 

dBM C ul.B^ise.J (3.1) 
(such a k exists by compactness). Note that 

k{s,r) = k{s/r,l), < r < s. (3.2) 

In other words, k{s, r) depends on (s, r) only through the ratio of s to r. 

Let X G M'^ and < r < s. Setting k = k{s,r), let ji, . . . ,jk be the first sequence of 
positive integers, according to the lexicographic ordering, such that (13.11) holds. That is, 
for < i < k — 1, having defined ji, . . . ,ji let jj+i be the first positive integer j such that 
{ Ji • • • ) Ji ) J } can be extended to a set of cardinality in N which satisfies (13. ip . 

For any x G M'^, if we set yi = x + scj-, then the points yi, . . . ,yk in dBs{x) satisfy 

dBs{x) c y^tlB:{y^)■ 

With yi, . . . ,yk defined thus, define 

F{x,s,r) := nl,{^ G N : y.(S°(y,) \ S,(x)) > 2}. (3.3) 

In the case x = we write simply F{s,r) for F{0,s,r). We think of F{x,s,r) as a set 
of configurations if containing a fence in the annulus Bs+r{x) \ Bs{x). We formalize the 
fence property as follows: 

Lemma 3.1. Suppose < r < s < 00 with s > 2r, and let x G M.'^. Suppose (p G F{x,s,r). 
Then for any z G ip\B°{x) we have Bp(^z,ip){.z) fl Bs-2r{x) = 0, and for any y G ipHBg^x) 
we have Bp^y^^){y) C B°^^,^.{x). 

Proof. Let z G Lp\B°{x). The line segment from z to x includes a point w in dBs{x). 
Then w lies in at least one of the balls B°{yi), and since we assume ip G F{x, s,r), there 
exists u E p n B°{yi) \ {z} such that \u — w\ < 2r. Hence by the triangle inequality, 
I2; — m| < I2; — w| + 2r. Therefore, since p{z,p) < \z — u\, the grain Bp(^z,<p){z) does not 
intersect Bs-2r{x), as asserted. 
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Let ?/ G (/9 n Bs{x). Take v in dBs{x), such that y lies in the hne segment from v to x. 
Since if G s, r), there exists t & ip \ {y} such that |t — f | < 2r. Hence by the triangle 
inequahty, piy,^) < \y — t\ < \y — v\ +2r, so that Bp(^y^^){y) C B°j^2r{^)i asserted. □ 

Recall that C M*^ is convex and compact with Vd(Vr) = 1, containing an open 
neighborhood of the origin, and Wn '■= n^^'^Wn for n G N. For proving results on point 
processes in Wn-, we introduce some further notation. Given < r < s, given G N and 
X G Wni we shall define -F„(x, s, r) similarly to s, r) but now with the fence involving 
only regions intersecting Wn- First we define the points yi = x + sej^ for 1 < i < k{s,r/2), 
similarly to the points y^ in the definition of F{x, s, r), such that 

dB.,{x) C Ufir^^)i?;/2(y.). (3.4) 

Let zi, . . . , Zk' be those y^ such that i?r/2(l/i) H 7^ 0. Note that k' < k{s, r/2). Set 

Fn{x,s,r) := nt,{^ G N : ^{B^izi)) > 2} n {v^ G N : ^ C (3.5) 

Lemma 3.2. Suppose y9 G N and G N (possibly with ip = (f). Let r,s,t > with 
s >2r and t > s + 4r, and x . Suppose y E ip H Bs{x) and z E ip\ Bt{x) Then: 

(i) If if G F{x^ s, r) and ip E F{x, t, r) then \y — z\ > p{y, ip) + p(z, ip). 

(a) If if E Fn{x, s, r) and ip G -F„(x, t, r) then \y — z\ > p{y, ip) + p{z, 

Proof. To prove (i), suppose ip G F(x,s,r) and if) G F{x,t,r). Then by Lemma 
13. H since t — s > 4r we have Bp{^y^^) C i?°^2r(^) Bpi^z,-^) H Bt-2r{x) = 0, so that 
Bp{y,^){y) nBp(^z^^){z) = 0, estabhshing part (i). 

To prove (ii), suppose instead that ip G Fn{x,s,r) and G F„(x,t, r). Let m and f 
be the points on the line segment yz such that — x| = s and |f — x| = t. Since Wn 
is convex, u and f are in Wn- Since m G dBs{x), by (13. 4p in the definition of F„(x, s,r) 
we have u G -Br/2(l/i) for some but then since also u G 1V„ we have y^ = Zj for some 
j < k'. Hence since ip G -F„(x, s,r), by (13.51) . there exists w E f H B°{zj) \ {y}, and by 
the triangle inequality 

\w — y\ < \u — y \ + \zj — u\ + \w — Zj\ < \u — y\ + 3r/2. 

Similarly, since v E Wn H dBt{x) and ip E Fn{x, t, r) we can find w' E ip \ {z} such that 
|w' — 2;| < \v — z\ + 3r /2. Hence by the hard-core property, 

p{y, ip) + p{z, ip) <\w — y\ + \w' — z\ < \u — y\ + \v — z\+ 3r, 

whereas since y,u,v, z are coUinear and |u — f| >t — s>4r, we have 

ly — z| > \u — y \ + \v — z\ + 4r, 

and part (ii) follows. □ 

We now give some probability estimates for the point sets F{x, s, r) and F„(a;, s, r). 

Lemma 3.3. There exists C3 G (0, 00) such that if < r < 2r < s and x E M.'^, then: 

(i) P(<l> ^ F(x, s,r)) < A;(s,r)c3exp(-cjV'^). 

(ii) If m,n E N and n > A and m E [n/2,3n/2] and x E Wn, then P(Xn,m ^ 
Fn{x,s,r)) < k{s,r)c3exp{—c^^r'^). 
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Proof, (i) With yi, . . . ,1/^ as in the definition f l3.3p of F{x, s,r), note that Vd{B°{yi) \ 
Bs{x)) > {bd/2)r^, so by subadditivity of measure, 

P($ ^ F{x, s, r)) < k{s, r)(l + (W2)r'^) exp(-(6V2)r'^), (3.6) 

and part (i) follows. 

For part (ii), first we claim that there is a constant c > 0, independent of n and x, 
such that if n > 4 and < r < n^/'^diam(Vr), then with k' and as defined just before 
(13. 5 p we have 

VdiBriz^) n > c"V^ l<t<k'. (3.7) 

Indeed, given i < k', if we choose y G Br/2{zi) fl W^n, then Br/2{y) C B.,.{z.i) so the claim 
follows from fICT]) . 

By (13. 7p . subadditivity and the fact that the binomial distribution is stochastically 
increasing in the success probability, we have 



m— ll 



nXn,m i i^n(x, r)) < k{s, r)[(l - rV(nc))'" + m(r7(nc))(l - r7(nc)) 

Hence by the inequality 1 — t < e~*, there is a further constant c' such that for all n > 4 
and m G [?t-/2, 3n/2], (so that in particular m — 1 > n/ A and m < 2n), for all x G and 
< r < s with r < 'n,-'^/'^diam(IV), 

nXn,m i i^n(x,s,r)] < fc(s, r) (1 + mr7(cn)) (1 - r7(cr2))"^-i 

< A;(s,r)(l + (2/c)r'^) exp(-(4c)"V'^) 

< A;(s,r)c'exp(-(c')"V). (3.8) 

Moreover, if r > n^/'^diam(iy) then s — r > n^/'^diam(iy) so k' = 0, and then trivially 
(13. 8 p still holds. This gives us part (ii). □ 

For G N and r > 0, define the set 

Sriv) ■=^xe^nBM0)\B2r(0)Six,(p). (3.9) 

Lemma 3.4. Let r > 0. Suppose (p G F{2r,r/2) fl F(4r, r/2) fl F[7r,r/2) , and suppose 
S*{ip) is bounded. Let (p^^ G N, (y^out G N 6e swc/i t/iat 

<^in C 52,(0) \ S:{p); <^o,t C R'' \ (B^riO) U 5;(^)). 

Then 

p{x, pU Pi,,) = p{x,Lp U Lfin^ Vont), X G (/^in U (v? H ^7, (0) ) ; (3.10) 
p(x, p> U v?out) = p(a;, pUpinU Pont), X G V^out ^ \ ^2,(0)). (3.11) 



Proof. First note that by Lemma [2.21 

p(x,v3Uv9in) = p(x,v2Uv9out) = p(a;, V^), X G (/? fl ^7, (0) \ 52, (0) . (3.12) 
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Hence, we can and do consistently define p'{x) for all x E ipU ipin U ipout, by 



p'{x) = p{x,LpU V^in) 

p\x) = p{x,ipU V?out): 



X e (^inU (^n 57,(0)); 
X e y^outU (¥?\52,.(0)). 



(3.13) 
(3.14) 



Assign grain radius p'{x) to each x G (pU(pin^{pout- We assert that with grain radii assigned 
in this way, each x E (p U (fin U v?out has a smaller grain- neighbour in ipU ipij^ U v?out, 

To verify this assertion, first suppose x G i?4r.(0). Then by the defining properties of 
the lilypond model, there exists yEifU (fin such that y is a smaller grain-neighbour of x 
under the radii p{-,ip U v?m)- Moreover, since ip G F(4r, r/2) fl F(7r, r/2), ii y ^ -87,(0) 
then by Lemma [3l2| we would have \x — y\ > p{x, ip U ipin) + p{y, ip U v^m), contradicting 
the statement that y is a grain- neighbour of x. Therefore y G -87,(0), so by fl3.13p we have 
p'{y) = p{y, ip U ipin) (and likewise for x). Therefore y is also a smaller grain-neighbour of 
X using the radii p'(-) as asserted. 

Now suppose instead that x ^ -84,(0). Then there exists y E ip U v^out such that 
y is a smaller grain-neighbour of x under the radii p{-,ip U v^out)- Moreover, since ip G 
-F(4r, r/2) fl F{2r,r /2), if y E -B2r(0) then by Lemma 13.21 we would have \x — y\ > 
p{x, if U v^out) + p{y, f U V'out) contradicting the statement that y is a grain-neighbour of 
X. Therefore y ^ -B2r(0), so by f l3.14p we have p'{y) = p{y, ip U v^out) (and likewise for x). 
Therefore y is also a smaller grain- neighbour of x using the radii p'(-) as asserted. 

We shall show that the radii p'{x),x E ip U ipm U v^out, have the hard-core property. 
This will suffice to give p'{x) = p{x, ip U v^in U v^out) for all x G V9 U (pin U (/9out as required, 
because as already mentioned, the lilypond model is the unique set of radii satisfying the 
hard-core and smaller grain-neighbour properties 

Let x,y be distinct elements of U y^in U v^out, with \x\ < \y\. We consider separately 
the case with \y\ < 7r, the case with |x| > 2r, and the case with < 2r and \y\ > 7r. 
These three cases cover all possibilities. 

In the first case with \y\ < 7r, both x and y are in ipi^ U {ip H -67,(0)). By (]3.13p and 
the hard-core property of the lilypond model on ipU yjin, we have p'{x) + p'{y) < \x — y\. 

In the second case with |a;| > 2r, both x and y are in yjout U {ip\B2r{0)). By fl3.14p and 
the hard-core property of the lilypond model on (y9 U v^out, we have p'{x) + p'{y) < \x — y\. 

Now consider the third case with \x\ < 2r and \y\ > 7r. In this case we have x E 
ipin U (v? n -B2r(0)) and y E yJout U (y? \ -67,(0)), so by (13.130 and fl3.14p . the assumption 
that ip E F{2r,r/2) fl F{7r,r/2), and Lemma [3. 2[ we have that p'{x) + p'{y) < |x — y\. 

Hence the radii p'{x),x E ipU ipi^ U y^out have the hard-core property as required. □ 

We have a similar result to Lemma [3.41 in the case of point processes in Wn- 

Lemma 3.5. Let r > and z E . Suppose ip E F„(z, 2r, r/2) fl F„(z,4r, r/2) fl 
F„(z,7r,r/2), with ipiW^) > 2. Let ipi^ E N, V^out eN be such that 



y^in C WnnB2r{x)\S;{^)- 



^ontCWn\{B-rr{x)US:{^)). 



Then 



p{x, ip U ipin) = P{X, V? U V^in U V^out), 
p{x, ip U (y9out) = Pix, ipUipinU V9out), 



a; G v^in U (v? n 57,(z)); 

X E (pont U (V5 \ B2riz)). 
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Proof. The proof is just the same as for Lemma 13.41 only using part (ii) instead of 
part (i) of Lemma 13.21 □ 

We now define U to be the set of all G N such that every point of (f has a unique 
smaller grain neighbour under the lilypond model based on ip, or y9(]R'^) < 2. 

Recall the definition ( j2l3l) of R{ip), G N. For a; G M'', and r > 0, n G N define 
subsets Er{x), Ur{x), Gr{x) and Gn,r{x) of N by 

Er{x) := {(^ G N : R{-y + ifi) < r, y E if H Bsrix)}; 

Urix) := {if EN : LP n B^r{x) E U}, 

and 

Gr{x) := Er{x)nUr{x)nn^j^^F{x,jr,r/2y, (3.15) 

GnAx) ■= Er{x)nUrix)nn%^Er,{x,jr,r/2). (3.16) 

For (y9 G N, we now define our radius of external stabilization i?ex(v^) by 

ReA^) := 9min{r G N : ^ G G^(0)}, (3.17) 

with min(0) taken to be +oo. The next result, in which we write Br for Br{0), shows 
that Rex has the external stabilization property. 

Lemma 3.6. Suppose r > 0. Then (i) if f E G^, then for any ip E N with iplB^r) = 0, 
we have: 



p{y,{^nBsr)U^) = p{y,{^''nB8r)Uij), y G U ((^ n Eg. \ 52r); (3.18) 

p{x,{(pnBsr)l^^) = pix,(pnBsr), X E (f H Bjr] (3.19) 

p{x, ((/?° n ^8^) u ^) = p{x, ip^ n ^8^,), X G n Bjr. (3.20) 

Also, (ii) if n eN, and x E Wn, and (p E G„^,.(x), then 

p{y,^) = piy,^""), yE^\B2r{x) (3.21) 



Proof, (i) Suppose ip E Gr- Then ip E F(2r,r/2) n F(4r,r/2) fl F(7r,r/2), and 
moreover ip E Er{0) so the set S*{ip) defined by (13. 9p is contained in B^r \ B,,.. Also 
S*{ip) = S*{ipnBsr{0)) by Lemma [2TT1 Therefore we can apply Lemma [3^ to ipHBgriO)] 
taking yjout = ^ and (pi^ = {0}, we obtain (13.181) from (13.111) and (I3.20p from (13.101) . and 
taking ipi^ = we obtain (13.191) from (13.101) . completing the proof of part (i). 

Part (ii) is proved by a similar argument, using Lemma instead of Lemma [231 □ 

The next result gives us tail bounds on _Rox(^)- 
Lemma 3.7. There is a constant E (0, oo) such that for all r > 0, 

P(<1> i a(0)) < C4exp(-C4 V'^/('^+i)) (3.22) 
and for n > 4 and m E [n/2, 3?2/2] and y E Wn, 

nXn,m i GnAv)) < C4 exp(-C4 V^/^'^+i)) . (3.23) 
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Proof. By Lemma 5.1 of [3], for all r > we have P($ G Ur{0)) = 1. Also, by the 
Palm-Mecke equation, see e.g. Theorem 9.22 of [TT], and f l2.4p . 

P($^E^(0))<E l{R{-x + ^)>r} 

= 6d(8r)'^P(i?($) > r) < bd{8r)'^ci exp(-crV'^/('^+^)). (3.24) 
Also, by Lemma 1X51 (i) and (13. 2p . 

8 

P(<l> ^ nf=iF(jr,r/2)) < A;(2j, l)c3 exp(-C3-V'^) 

i=i 

and combined with (I3.24p this gives (I3.22p for a suitable choice of C4. 

Next, suppose n > 4 and n/2 < m < 3n/2, and y G VFn- Then by Lemma [2.5[ 

P(Xn,m ^ ^r(l/)) < " / P(i?(-X + Xn,m-l) > r)dx 

< (3/2)c2exp(-c^V'^/('^+i)). (3.25) 
Combined with Lemma 13.31 (ii) and (13. 2p , this gives us (I3.23p . □ 

4 Sub-exponential decay 

In this section we give tail bounds on the distribution of the size of the component of the 
Poisson lilypond model containing a typical Poisson point, with 'size' measured either by 
cardinality, or by metric diameter, or by volume. 

For G N, define as follows the directed graph Q{(p) = {(p,E{{p)) with vertex set (p 
and edge set E{ip). A pair {x,y) is in E{ip) if ?/ is a smaller grain-neighbour of x. Let 
^*((/?) denote the associated undirected graph. For x G v^, let C{x,ip) denote the cluster 
at X, that is, the set of points of ip that are connected to x by a path in the undirected 
graph Q*{^p). Let C"(x, ip) denote the union of lilypond grains centred at points of C(x, v?), 
i.e., the connected component containing x of the set [Jx^ipBp(^x,ip){.x). 

Theorem 4.1. There are strictly positive constants C5, ce, cy such that 

C5^exp(-C5r'^) < P(diam(C"(0,$O)) > r) < C5 exp(-C5 V'^/('^+^)), r > 0; (4.1) 
C6-^exp(-C6t) < P(Kz(C"(0,$0)) >t) < C6exp(-C6Hi/('^+i)), t > 0; (4.2) 
cf^exp(-C7n2) < P(cardC(0,$°) > n) < C7exp(-cf ^n'^/('^+^)), n > 0. (4.3) 

Note that in each of (14. ip . (14. 2 p and (14. 3 p the power of r, t or n in the exponent is 
different in the lower bound than in the upper bound. It is an open problem to make 
these bounds sharper. 

Sharper bounds are available in the analogous setting for lattice and continuum per- 
colation. Consider for example the geometric graph on with each pair of points 
connected by an edge if and only if the distance between them is less than a constant r*, 
with r* chosen to be subcritical. Then results like (14. 3 p and (14. ip hold with exponents of 
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the form c-n, respectively c-r, in both the upper and lower bound, though not necessarily 
both with the same c (see Section 10.1 of [H]). A result like fl4.ip holds for subcritical 
lattice percolation with exponents of the form c-r; see (6.10) of [Tj. This implies a bound 
like ( 14. 2 p for the Boolean model associated with the subcritical geometric graph just men- 
tioned, with c ■ t in the exponent. A similar lower bound also holds. 

The next lemma will be used in proving Theorem 14.11 and again later. Recall the 
definitions f l3.15p and f l3.16p of Gr{x) and Gn,r{x) respectively. 

Lemma 4.2. Let r > and x G W^. 

(i) If (f E Gr{x), then C{y, ip) C B^r{.x) for all y E B^rix), and 



(a) If n E N, and x E Wn, and if E Gn,r{x), then C{y,(p) C i?5,.(x). 

Proof, (i) First, we assert that each y E ip (1 B^r{x) has a unique smaller grain 
neighbour in ip. Indeed, y does not have any grain- neighbour (in either (p or ipd Bgr{x)) 
outside Bjrix), by Lemma [121 because ip E F{x,5r,r /2) fl F(x, 7r, r/2). Also, p{u,ip) = 
p{u,{p n Bgr{Qx)) for all u G n Bsr{x) by Lemmas 12.11 and 12. 2^ because ip E Er{x). 
Hence, the unique smaller grain neighbour of y in (p (1 Bgr{x) (which it has because 
(p E Gr{x) C Ur{x)) is also its unique smaller grain neighbour in (p, justifying the assertion. 

Hence, in the graph Q{ip), each vertex inside B^rix) has an out-degree of 1. This 
implies that for any path in Q*{'P>) starting inside B^r{x) and ending outside B^r{x) but 
with all vertices except the last inside B^r{,x), if the direction of the edges in Qi^p)) is 
taken into consideration the path can reverse its direction at most once. That is, such a 
path must consist of a directed path in the forward direction (possibly of zero length), 
followed by a directed path in the reverse direction (also possibly of zero length). Hence, 
if ?/ G n B^j.{.x) and C(?/, ip) is not contained in i?5r(x), then either there is a descending 
grain-chain in ip starting inside B2,r{x) and ending outside i?4r(x), or there is a descending 
grain-chain in tp starting outside B^r{x) and ending inside i?4r-(a;). 

Since G Er{x), we have no descending chain starting inside B^r{x) and ending 
outside Bir{x) or starting in B^ri^x) \ B^rix) and ending inside B^rix). Moreover, since 
ip E F{x, 4r, r/2) nF(x, 6r, r/2), by Lemma [3l2] there is no edge of G{ip) with one endpoint 
outside BQr{x) and the other endpoint inside B^rix). This shows that C{y,(p>) C B^r{x), 
and since also ip E F{x,5r,r /2), by Lemma [3.11 we have (14. 4p . 

The proof of (ii) is similar. □ 

Proof of l \4-ll ■ By Lemmas 14.21 and 13. 7[ we have for r > that 




(4.4) 



P({C"(0,$' 



) C BUO)r) < P($ ^ ^,(0)) < C4exp(-c; 



-1 d/(d+i) 



(4.5) 



-4 



and the upper bound in fl4.ip follows. 

For the lower bound, observe that for all r > 0, we have 



P[diam(C"(0, $°)) > r] > P[$(52,(0)) = 0] = exp(-6d(2r)'^). □ 
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Proof of i \4-3^ ■ Choosing r so that 2bd{Qr)^ = n, by (14. 5 p and a standard Chernoff-type 
tail estimate for the Poisson distribution (see e.g. Lemma 1.2 of [TT]) there are constants 
c, c' such that for n > 1 we have 

P[card(C(0,$°)) >n]< cexp(-c-V<^/('^+^)) + P[<1>°(56,(0)) > 2bd{6rY] 

<c'exp{-{l/c'y/^''+^^) 

and the upper bound in (14.30 follows. 

Let e be the unit vector (1,0,..., 0) in W^. For i e N, let Bi be the closed ball in W^, 
centred on 9~*e and having radius 9~*~^. Since Q"*-*"^-' — 9~* = 8(9~*), for alH > 2, x G -Bj 
and y G -Bj-i we have 7(9~') < \x — y\ < 9"*^*""'^^ and also |x| < 10(9~*~^). 

Observe that there is a positive constant c such that for all n G N we have 

P({$(59(0)) = n}n rt^dHBi) = 1}) = (^{^9-^'^'^')^ exp{-b,9') 

> c^^ exp(— cra^). (4.6) 

If the event inside the left hand side of (14.61) occurs, then labelling the point of $ fl 
as Xi, each point Xi for 1 < i < n has a smaller grain-neighbour in to its left in the 
collection {0, xi, X2, ■ ■ ■ , Xn}- Indeed, every point in this collection to the right of Xi, and 
also every point of $ \ -89(0), is distant more than 2D{xi, $") from Xi, so cannot be its 
smaller grain neighbour. 

Therefore, if the event inside the left hand side of (14. 6 p occurs, then all 
in C(0, $"), so that in this case cardC(0, $") > n + 1. Hence, the lower bound in (14. 3 p 
holds by (jM]). □ 

Proof of 1 14-2^ ■ Taking r so that bar'^ = t we have by (14. ip that 

P[Ki(C"(0,$°)) >t]< C5exp(-c^V'^/('^+i)) = C5exp(-c^i(t/6d)i/('^+i)) 
while for the lower bound, we have 

P[Vd(C"(0, <1>°)) >t]> P[<1>(52,(0)) = 0] = exp(-2'^t). □ 



5 Central limit theorems 

In this section we derive central limit theorems (CLTs) associated with the lilypond model. 
We consider a sequence of binomial or Poisson point processes with finite total number 
of points over an expanding sequence of bounded regions in R''. Thus we consider, for 
example, Vd(Z($„)) as defined in Section^ rather than Vd{Z nWn) ■ Our choice means 
that we can directly apply results in [l^ or [13], although it should be possible to obtain 
similar results for Vrf(Z($) n Wn). 

For n,m G N, recall the definition of Xi, X2, . . . and Xn,m at (^M)- Set Wn := n^/'^W. 
We consider the restricted Poisson process $„ := $ fl Wn, and also the binomial point 
process Xn '■= Xn,n- For measurable g : M+ — ?■ M, we give CLTs for sums of the form 
Hg{^n) and Hg{xn), where for all (/? G N with 2 < (p(R^) < 00 we define 
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and if (^^(R'^) G {0, 1} we set Hg{ip) := 0. For example if g{t) = baf^, then by the hard-core 
property Hg{^n) is the total volume Vrf(Z($„))l{$(vt^^)>2} and Hg{xn) is Vd{Z{xn)) for 
n>2. 

Theorem 5.1. Suppose that there exists finite /3 > such that g : M.^ — t- M satisfies the 
growth bound 

snp {1 + t^)-^\g{t)\ < oo, (5.1) 

and that g is not Lebesgue- almost everywhere constant. Then there exist constants < 
Tg < < oo (dependent on g but independent of the choice of the convex set W) such 
that as n ^ CO we have n~^Var(_f/'c,($„)) — cr| and n^^'Vax{Hg{xn)) Tg, and 

n-^'\Hg{^r.) - ^Hg{^r.)) A iV(0, aj); (5.2) 
n-'l\Hg{xn) - EHgixn)) A N{0, t^). (5.3) 

The next theorem provides a similar CLT for the number of components. Recall that 
/t(A) (for A C M"^) denotes the number of connected components of A] for G N with 

2 < ifiW^) < oo we define 

which is also the number of components of the graph ^* (</?), and if y9(R'^) = 1 we set 
H^{^) := 1, and if ^{W^) = we set H^{^) := 0. 

Theorem 5.2. There are constants < < < oo such that n~^Var(i7«;($„)) — )■ o"^ 

and n^^Yai{H t^{xn)) '^'i, o,nd 

n-^/\H^{^n) - Ei/.($„)) A iV(0, al)- (5.4) 
n-^'\H^{Xn)-^H^{Xn)) -A iV(0,r2). (5.5) 

We shall prove Theorems 15.11 and 15.21 by using results from [14J which have the ad- 
vantage of showing that the limiting variance is non-zero. An alternative would be to use 
results from [13] (an approach based instead on sub-exponential stabihzation) which could 
be used to give a formula for cx^ in terms of integrated two-point correlation functions, 
and also to provide Gaussian limits for the random measures associated with the sums 
Hg{^n) and Hg{xn)', moreover, using Theorem 2.5 of |15], it should be possible to provide 
Berry- Esseen type error bounds associated with (15. 4p converging to zero at rate 0{n^^^^^) 
for any e > 0. However, we do not give details of these alternative approaches here. 

The approach of [14j is based on a notion of external stabilization. Given a real- valued 
functional H{if ), defined on finite G N in a translation-invariant manner (i.e. with 
H{x + if) = H{x) for all nonempty y9 G N and all x G M.'^), we say that a nonnegative 
random variable i? is a radius of external stabilization for H if there is a further random 
variable A such that 

i/(($° n 5^(0)) utp)- H{{^ n 5^(0)) u ^) = A 

for all ip E N with '?/'(M'^) < oo and ilj{Bj^{0)) = 0. If if has a radius of external 
stabilization R with P(i? < oo) = 1, we say H is externally stabilizing and refer to A as 
the add- one cost of H. 

To prove Theorems 15. II and 15.21 we shall use the following result, which is Theorem 2.1 
of [H]. 
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Theorem 5.3. Suppose H is externally stabilizing with add-one cost A, and satisfies the 
moments condition 

sup sup E{\H{xlJ - H{Xn,m)\') < oo (5.6) 

n,mGN:n>4,me[n/2,3n/2] a;GVK„ 

along with (for some [5 G (0, oo) ) the growth hound 

\H{v)\ < /3(diam(v9) + </^(M^))^, ¥^ e N, < ^^(M'^) < oo. (5.7) 

Then there are constants < < cr^ < oo such that n~^Var(if ($„)) — )■ cr^ and 
n~^Vai{H{xn)) — ^ t"^, and 

_EiJ($„)) A iV(0,a2); 

n-i/2(/7(x„)-Ei7(x„)) AiV(0,r2). 

Moreover, if A /ias a non- degenerate distribution then > 

Proof of Theorem \5.1[ We shall use Theorem 15.31 Recall the definition 13. 171 of -Rex(v^)- 
For G N with Rex{'^) = 9r < oo, set 

By Lemma [3^ if G N with -Rex(v') < C)0 and G N with ?/'(M"') < oo and 'ipiBR^^(^){0)) = 
0, then 

i7,((v9° n i?«..(^)(0)) U ^) - if,(^ n Bn^^i^m) = A,((p). 



Hence, since i?ex($) < oo almost surely by Lemma [321 Hg is externally stabilizing. Using 
(15. ip . it is easy to see that H = Hg also satisfies the growth bound (15.71) . We need to 
check that Hg satisfies the moments condition (15. 6p . and that the add-one cost Ag($) is 
non-degenerate. We demonstrate these in Lemmas 15.41 and 15.51 below. Given these, we 
can apply Theorem 15.31 to get the result. □ 



Lemma 5.4. Suppose g : M_|_ — )■ M satisfies the growth bound (15. iP . Then the functional 
H = Hg satisfies the moments condition (\5.6h . 

Proof. Since Xn,m C. Wn, by Lemma \3M (ii), and Lemma [3?7l for all r > 0, n > 4, 

X G Wn and m G [n/2, 3n/2] we have 



P 



Hg{Xn,J - Hg{Xn,m) 7^ \ M^, \ ^ 9{p{y,f)) 

Z&Xn.m<^B2r(x) j yeXn,rnnB2r {x) 

< C4exp(-C4V'^/(^+i)). (5. 



Now observe that there is a constant Cg such that by the assumed growth bound (15.11) 
on g and the bound p(x, (p) < D{x, 99), if 2 < {p{B2r{x)) < 4:bd{2rY, then 



Y 9{p{z,^'')) 1 - Y 9{p{y.v)) 



< C8(l + r 



(5.9) 
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By fl2.1UI) . if n > 4 and m G [n/2, 3n/2] and x e Wn and < r < n^/'^diam(W^), then 

Exn,m(52.(x)) = im/n)VdiB2r{x) n G [i2co)-\2rY, 2bd{2r)''] . (5.10) 

In particular, if {ScqY^'^ < r < r;,^/'^diam(iy), then ¥.Xn,m{B2r{x)) > 2. By a standard 
Chernoff-type tail estimate for the binomial distribution (see e.g. Lemma 1.1 of pjj), and 
f IS.lOp . there is a constant Cg G (0, oo) such that if (Scq)^/'^ <r< n^/'^diam(VI^), then 

1 - F[2 < Xn,m{B2r{x)) < 46,(2r)'^] < 2exp(-C9-V'^). (5.11) 

Moreover, if r > n^/'^diam(W^) then since Vd(H^) = 1, the Bieberbach (isodiametric) 
inequality (see e.g. [llj) yields 

^hd{2rY > 2'^'^+'^nbd(diam{W)/2y > 2^'^+^n, 

so if 2 < m < 3n/2 then trivially f l5.1ip holds in this case too. Combining (15. 8p . (15. 9 p and 
(15.1 ip gives us, for all n > 4, m G [n/2, 3n/2], x G Wn and r > {ScqY^^, the tail bound 

P[l^.(x:,„.) - Hg{xn,m)\ > cs{l + r'+^)] < C4exp(-C4-V'^/('^+i)) + 2exp(-C9/), 
which suffices to give us the moments bound (15. 6 p for Hg. □ 

Lemma 5.5. The random variable Ag($) has a nondegenerate distribution. 

Proof. By assumption g is not almost everywhere constant, so by the Lebesgue density 
theorem (as stated in e.g. |TT]) there exist numbers a > 6 > and ti > 0, ^2 > and 
e G (0,min(ti,t2)/10) such that f^^^' l{g{t) > a}dt > 0, and J^^~^' l{g{t) < b}dt > 0. 
Choose r > max(ti,t2) such that P($ G F(99r, r)) > (this is possible by (l3l6l)). 

Enumerate the points of $ as Xi, X2, . . . with < IX2I < ■ ■ ■ . Set D := \Xi — X2I 
and Rq := inf{|x| : x G Bd/2{Xi)}. For i = 1,2 define the event 

Ei := {\D\ <e}n {Ro G {ti,ti + e)} n {IX3I > 99r} n {$ G F(99r,r)} 

Then W{Ei) > and F{E2) > 0. Also, given Ei occurs the value of p(0, $'^) is equal 
to Rq and Hg{^^ n 5^(0)) - Hg{^ n 5^(0)) = g{Ro) for any s > lOOr, so (provided also 
i?ex($) < 00) the value of Ac,($) is equal to g{Ro). Moreover, given E^, the distribution 
of Rq is absolutely continuous on (tj,ti + e) with a strictly positive density. Therefore 
P[A<,($) > a\Ei] andP[A<,(<l>) < b\E2] are both strictly positive. Hence P(A3(<I>) > a) > 
and P(Ag($) < 6) > 0. Thus Ac,($) has a nondegenerate distribution as required. □ 

We now proceed towards a proof of Theorem 15.21 We wish to show the the functional 
has an almost surely finite radius of external stabilization. In fact we shall again use 
i?ex5 as defined at (I3.17p . 

Suppose -Rex(v') = 9r < 00. Let ip* := ^xeipnB3r{0)C{x,(p). Let N{ip) be the number of 
components of the subgraph of induced by the set of vertices ip*. Let No{ip) be the 
number of components of the subgraph of ^(v2°) induced by ip* U {0}. Let 

A,(<^) :=iVo(<^)-iV((^). 

The following lemma says that -Rex serves as a radius of external stabilization for the 
functional H^: 



18 



Lemma 5.6. (i) // ^9 G N with -Rex(v^) < oo and E 'N with iIj{M.'^) < oo and 
V'(5fi,,(^)(0)) =0, then 

H,{{^' n i?K..{^)(0)) u ^) - H^iiv n i?ij,,(^)(0)) u v^) = A.(^). (5.12) 

(a) Let n eN. Suppose x G Wn, and if G Gn;r{x). Then 

Proof. Suppose -Rex(v^) = 9r < cxd. Then by Lemma [3.6^ the radii of hlypond grains 
centred outside i?2r(0) are unchanged when a point is inserted at 0, while the radii of 
grains centred inside B-jriS^) are unaffected by changes to ip outside Bi^riS^). 

Set ifi := {(p n Bgr{0)) U ip. Let x E ip* and y E (fi \ (f*. Then we claim that 
{x,y} ^ and {x,?/} ^ Q*{ip^). Indeed, by Lemma W72\ we have x G 55^(0), so for 

y ^ i?7r(0), the claim follows from Lemma and the fact that ip G F{7r, r/2)nF(5r, r/2) 
(so also G F(7r, r/2) fl F(5r, r/2)). In the case ?/ G i?7r(0), the claim follows from 
the definition of ip*, along with the fact that p{y,ipi) = p{y,(p) and p{x,ipi) = p{x,ip). 
Similarly, {0,y} ^ g*{^',). 

Therefore the components of Q*{ipi) induced by all y G ipi \ ip* do not meet the 
components containing all x E ip* and this still holds after the addition of 0. So the 
contribution of these components to H^{(p\) — Hf^^ipi) is zero and (15.121) follows. 

Part (ii) is proved by much the same argument as for Part (i), using Lemma [375] instead 
of Lemma 13.41 and Lemma 13.21 (ii) instead of Lemma 13.21 (i). □ 

Proof of Theorem By Lemmas 13.71 and 15.61 is externally stabilizing. It is 
easy to see that if^ satisfies the growth bound (15. 7p with [3 = 1. We need to check that 
i^re satisfies the moments condition (15. 6p . and that the limiting add-one cost Ak($) is 
non-degenerate. We demonstrate these in Lemmas 15.71 and 15.81 below. Given these, we 
can apply Theorem 15.31 to get the result. □ 

Lemma 5.7. The functional H = satisfies the moments condition (15.61) . 

Proof. By Lemmas 15.61 (ii) and 13. 7[ for n, m G N with n > A and n/2 < m < 3n/2, 
and X G Wn and r > 0, 

Therefore, for all £ G N, 

ni^.(x:,J-^.(Xn,JI >^) <cexp(-c-V/('^+i))+P(x:,„(i?7r(x)) >£). (5.13) 

Now take r := r{£) := (l/7)(£/36d)^/'^. Then ExnABjr) < fed(3/2)(7r)°' = £/2. Hence the 
last probability in (I5.13P decays exponentially in £ by standard Chernoff tail estimates for 
the binomial distribution (see e.g. Lemma 1.1 of [H]), so that overall there is a constant 
c such that the right hand side of ( 15.13P is bounded by cexp(— c"^^^''*^'^"'"^-'), for all £ < n, 
independently of n and m. Moreover, the left side of (I5.13P is zero for £ > n. This gives 
the uniform bound (15. 6p on fourth moments for H^- 
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Lemma 5.8. The limiting add-one cost Ak($) has a non- degenerate distribution. 

Proof. Let e G M*^ be a unit vector. Enumerate the points of $ as Xi,X2, ■ ■ . with 
|Xi| < IX2I < ■ ■ ■ . Given r > 1, define events i?o,r and Ei^r as follows: 

Eo,r := {$(5i(5e)) = 2} n {IX3I > 99r} n {$ G F(99r,r)}; 

Ei^r ■■= {$(5i(0)) = 1} n {<l>(5i(5e)) = 2} n {IX4I > 99r} n {$ G F(99r, r)}. 

Then by (13.61) . for large enough r we have F{EQ^r) > and P(i?i_r) > 0. But if Eq ^ occurs 
and -Rex('^') < 00 then A^{^) = 0, while if Ei^r occurs and Re^i^) < 00 then A^{^) = 1. 
This gives us the result. □ 



6 Percolation theory for the enhanced model 

One interpretation (in = 2) of the absence of percolation in the Poisson lilypond model, 
is that a frog is unable to travel infinitely far by a continuous path along the lily pads. 
For the sake of greater realism, it is natural to ask what happens if the frog is allowed to 
jump. More mathematically, for 5 > we consider the enhanced union set Z^, where for 
any A C M'^ we set := UxeABs{x) so that in particular 

Z' := U i?,(,,$)+,(x). (6.1) 

We investigate the connectivity properties of Z^ (the parameter S represents half the dis- 
tance which the frog is able to jump). In particular, we are concerned with the probability 
that Z^ has an unbounded connected component. Given S, the event that this occurs is 
invariant under translations of the Poisson process $, so by the ergodic property of this 
Poisson process (see Proposition 2.6 of [lO]), or alternatively by an argument using the 
Kolmogorov zero-one law, this probability is either zero or 1, and it is one if and only if 

P(Eoo(5)) > 

where E^dS) denotes the event that there is an infinite component of Z^ containing the 
origin. Accordingly, we define the critical enhancement 6c '■= Sc{d) by 

4 := inf{5 > : P(Soo(5)) > 0}. 

Our main result in this section says that if the range of our jumping frog is sufficiently 
small then it is still unable to travel infinitely far. 

Theorem 6.1. If d = I then Sc = 00. If d > 2, then < Sc < 00. 

For d = 1, it is easy to see that 6c = 00. For d > 2, the fact that 6c is finite is 
immediate from the basic fact in continuum percolation, that for sufficiently large r the 
union of balls Uxe'S?Br{x) percolates (see for example [7] or [ID]). Thus, to prove Theorem 
16.11 we need only to consider the case with d>2 and show 6c > 0. 

Let Dr{x, (p) be the minimum of all non-zero pairwise distances between lilypond balls 
centred in Bsr{x). That is, let 

Dr{x, ip) = min{|M — v\ — p{u, ip) — p{v, p) : u,v E p H Bsr{x), \u — v\ > p{u, p) + p(f , p)} 

and set Dr{x,p) = +00 if there are no such pairs {u,v). Note that Dr{x,p) is strictly 
positive since, if finite, it is the minimum of a finite set of strictly positive numbers. 
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Lemma 6.2. Given x G M'^ and r > 0, the event 

{$ G Gr{x)}r}{Dr{x,^) > 6} 
is measurable with respect to o"($ fl Bgr{x)). 

Proof. Let G N. By Lemma [271] S{y,Lp) C -8, (2/) if and only if S{y,ip fl Br{y)) C 
Br{y). Hence, ip G (^^(a;) if and only if fl -B9r(a;) G (^^(x). If this is the case, then by 
Lemmas 12.11 and 12.21 we also have Dr{x,ip) = Dr{x,ip fl Bgr{x)). Hence the displayed 
event is identical to the event 

{($ n Bgrix)) G Grix)} D {Drix, $ H Bgrix)) > 6} 

which is measurable with respect to (t($ fl Bgr{x)). □ 

A key claim is the following: 

Lemma 6.3. Let 6 G (0, 1/2) and let r > 2. Let x G M.'^. If there is a continuous path in 
from R'^ \ B^rix) to Br{x), then the event {$ ^ Gr{x)} U {Dr{x, $) < 26} occurs. 

Proof. Suppose that {$ ^ Gr{x)} U {Dr{x,^) < 26} does not occur. Let T denote 
the union of components of Z which intersect with B2r{x). Any component of Z with 
all its Poisson points outside B^rix) is contained in M.'^ \ B2r{x), by Lemma [3TT] because 
$ G F{x,3r,r/2). Hence by Lemma [4. 2 [ 

T C Uye^nB,A-)C'iy, C BM. (6.2) 

Consider all lilypond balls centred at Poisson points outside T. Those centred inside 
Bsr{x) are at distance more than 26 from T, because Dr{x,^) > 26 by assumption. 
Those centred outside Bsr{x) do not intersect Bjr{x) because of the assumption that 
$ G F{x,8r,r/2) and Lemma |3.H and so are also distant more than 26 from T, since 
flOj) holds and r > 26. 

Thus, the set T is distant more than 26 from Z\T, and hence is disconnected from 
the rest of Z^. Finally, by definition of T, C Bjr{x) and Z^\T^ is disjoint from Br{x), 
both because 5 < r, so that there is no continuous path in Z^ from Br{x) to M'^ \ Bjr{x). 
□ 

Proof of Theorem \6.1[ Recall the definitions of R{'^) and Gr{x) at (12. 3p and (I3.15P 
respectively. By Lemma ISTl P($ ^ G'r(O)) — t- as r — )■ oo. Also, Dr{x) is a strictly 
positive random variable. Hence, given e > 0, we can choose r > 2 large enough and 
5 G (0, 1/2) small enough such that 

^Gr{0)}U{Dr{0) <26}) <e. (6.3) 

Note that P({$ ^ Gr{x)} U {L'r-(x) < 25}) is the same for all x. 

Now divide M*^ into boxes (cubes) of side 2rd~^^'^, labelled Qz,z G Z'^, by setting 

= 2rrf-i/2^ + [-rd~^/^rd~^/Y- 



21 



Define the random field {Y^, z G Z'^) by 

Yz ■= i — l{$^Gr(2rd-V22)}u{Dr{2rd-i/22,$)<5}- 

Since Qz C Br{2rd~^^'^z), by Lemma [6.31 if tliere is a continuous patli in from to 
M*^ \ Ber{2rd-'^/'^z), then = almost surely. 

If there is an infinite component in Z^, then there must be an unbounded continuous 
path in Z^, and by taking successive boxes along the path, there is an infinite sequence 
{zi, Z2, -23, . . .) of distinct elements with each Zi G Z'^ and \\zi — 2;i+i||oo < 1 for each i, such 
that Yz = for all i. 

Given z G Z*^, by Lemma 16.21 the random variable Y^ is measurable with respect to 
o"($ n i?9.r(2r(i~^/^2)). Hence the random field {Y^^z G Z"^) is independent of Y^i for all 
all sites z! with 2rd~^l'^\z' — z\> 18r, i.e. with \z' — z\> ^d}!'^ . Thus Y;, is independent 
of Y^i whenever the graph distance between z and z' exceeds 9(i. In fact, (F^, z G Z"') is a 
9(i-dependent random field in the sense of [7]. 

Let j9c be the critical probability for site percolation on the lattice with vertex set Z'^ 
and edges between each pair [z, z') with \\z — z'\\oo = 1; it is well known [7] that pc > 0. 
By ([7], Theorem (7.65)), and (16. 3p . we can choose r to be so large and 6 to be so small 
that the random field (Yz, z G Z"^) stochastically dominates a random field {Y^, z G Z'^) 
consisting of independent Bernoulli random variables with P(K,' = 0) = Pc/2 for each z. 
Thus, with this choice of 6 there is almost surely no infinite path through the lattice of 
sites with Yz = 0, and hence no infinite component in Z^; hence (5c > as asserted. □ 
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